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© Micro- and nano-porous metallic structures. 



© Micro- and nano-porous metallic structures are produced by contacting an alloy of at least two metals having 
different electrochemical activity in predetermined proportions with an electrolyte and applying an electrochemi- 
cal potential to said electrolyte wherein said alloy acts as the anode until all or a portion of the more 
electrochemically active metal is dissolved from the alioy leaving porous metal. The resultant structure has 
substantially uniform and interconnected porosity in pore sizes ranging from 3 nanometers to 10 micrometers. 
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MICRO- AND NANO-POROUS METALLIC STRUCTURES 

The present invention relates to certain porous metallic structures and processes for making same. The 
invention particularly relates to metallic structures having very small pore sizes being of uniform and 
interconnected porosity. The porous metallic structures have utility, for example, as metal membranes or 
filters. 

5 Broadly described, the porous metallic structures are produced by mixing at least two metals having 

significantly different electrochemical activities, according to predetermined proportions, to form a multi- 
metallic alloy; contacting the alloy with an electrolyte; and applying an electrochemical potential between 
the metal and the electrolyte wherein the alloy acts as an anode until all or a portion of the more 
electrochemically active metal is dissolved from the alloy leaving porous metal. The porous metal in another 

id embodiment of the present invention is heat annealed i.e. coarsened to increase the pore size and 
mechanical stability of the porous metal. The structures are produced with uniform porosity in pore sizes 
ranging from about 3 nanometers to 10 micrometers. The porous structure can be produced in thicknesses 
ranging from 0.0001 to 10 mm. Depending upon the elemental composition of the membrane, the volume 
fraction of porosity attainable in the structures ranges from 30-80% and all of the porosity is substantially all 

75 interconnected. Membranes prepared in this way are mechanically tough, chemically inert, and heat 
resistant (e.g. thermally sterilizable). Such a membrane may be employed as a filter to separate particles, of 
size down to the molecular level, from gases and liquids. 

Further objects and advantages of the present invention will be better understood by reading the 
following detailed description of the presently preferred exemplary embodiments of this invention in 

20 conjunction with the accompanying drawings, of which: 

FIGURE 1 shows polarization curves of Ag-Au alloys in 1M perchloric acid, (a) pure Au (b) 
Ago.5DAuo.50 (c) Ago.ooAuo.4o (d) Ag 0 . 70 Au 0 ^o (e) Ago.74Auo.26 (f) Ag 0 .8oAu 0 . 2 o; 

FIGURE 2 shows the electrochemical cell configuration for membrane production; 

FIGURE 3 shows curves for the normalized pore size distribution for Ago.74Au 0 . 26 alloy annealed 10 

25 minutes at the indicated temperatures. The pore size area, A, was normalized using the maximum pore area 
which developed at each of the annealing temperatures, A max (temp " C). These parameters are, A^SOO)- 
= 2.45 urn 2 , A max (700) = 0.95 urn 2 , A max (600)=0.55 urn 2 , A max (500) = 0.28 urn 2 and A max (400) = 0.24 urn 2 ; 
and 

FIGURE 4 shows curves for normalized pore aspect ratio distribution for Ag 0 .74Au 0 ^ 6 alloy annealed 

30 10 minutes at the indicated temperatures. 

The porous metal structures are produced starting with an alloy of prescribed composition using a 
combination of electrochemical and thermal treatments. The starting composition of the alloy determines 
the volume fraction of porosity and the minimum pore size available from the particular alloy. The starting 
composition is determined using a scheme based upon percolation theory. The alloy is composed of 

35 elemental metals of prescribed composition with a large difference in electrochemical activity. The relative 
electrochemical activities of the elements are determined by the difference in the thermodynamic equilib- 
rium metal/metal ion electrode potentials, AE*. as given by the electromotive force series of standard 
reduction potentials. By the application of an appropriate electrochemical potential in a suitable electrolyte 
the reactive element is selectively dissolved out while leaving the more noble element in the form of a 

40 sponge or membrane retaining the initial shape and dimensions of the sample. This allows the shape of the 
membrane to be pre-machined prior to the corrosion and thermal treatments as discussed below. 
Depending upon the system, the dissolution process may take place at various temperatures. 

For a particular alloy system the "de-alloying threshold, p\" is defined based upon the difference in the 
equilibrium metal/metal ion electron potentials of the elements in the alloy. The p. defines the minimum 

45 possible concentration of the reactive element in the alloy necessary for forming the porous membrane 
from the remaining noble-metal and is determined from percolation theory in the following way. For 
elements in a binary alloy with a AE* greater than 0.3 Volt and less than 1.0 Volt, selective corrosion of the 
more reactive element can only proceed by a two-dimensional percolation process. This requires a 
minimum composition of the more reactive element of ca.60% atomic fraction (at%). For elements in a 

50 binary alloy with a AE° of 1.0 volt or greater, selective corrosion of the more reactive element proceeds by 
a three-dimension percolation process which allows for a minimum composition of the more reactive 
element of ca.30 at %. The atomic concentration of reactive element determines the volume fraction of 
porosity achievable in the membrane and the process of the present invention enables a volume fraction 
range of 30-80%. Thus, the metal to be removed comprises 30-80% of the binary alloy depending on the 
AE* of the alloy. The 80% upper limit is a result of the requirement that the remaining noble porous metal 
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be able to span the geometry or shape of the structure. For starting compositions greater than 80 at % of 
reactive element the remaining porous metal cannot maintain the initial shape of the structure and 
collapses. The electrolyte is chosen such that it has high solubility for the dissolved reactive metal 
component. The electrochemical potential is chosen to be equal to (or as close to, as allowed for by the 

5 electrolyte) the equilibrium metal/metal ion electrode potential of the more noble metal element comprising 
the alloy, assuming a metal ion concentration of ca.10~ e M. 

Metal removal by the corrosion procedure is in one embodiment of the invention, carried out to a 
prescribed depth as a fraction of the sample thickness, resulting in porosity to that depth. 

The porous metallic structure resulting from the corrosion process, may be reinforced by two additional 

70 procedures. The resulting porous structure may be reinforced by applying i.e. lacquering a tightly adhering 
material impermeable to the electrolyte to selected areas of the initial alloy before contacting the alloy with 
the electrolyte and initiating the corrosion procedure. This prevents subsequent dissolution of the reactive 
metal component (e.g. silver in a silver-gold system) from the covered areas. The undissolved metal i.e. 
original solid, ductile alloy, is present between the dissolved regions thus providing structural support to the 

75 membrane. The solid portions of the porous structure allow for clamping and support of the structure in 
applications. 

Another method for structurally supporting the porous metallic membrane resulting from the corrosion 
process, is to place a metal (e.g. copper) macroscopic mesh or other design on ail or a portion of the 
surface of the porous metallic structure. The reinforcing material forms a composite structure providing 
20 great structural integrity. 

In an embodiment of the present invention, after the corrosion procedure the resulting porous metallic 
structure may be subjected to heat annealing (coarsening) to increase the pore size, mechanical stability, 
and ductility of the structure. The unannealed structure is typically placed in ,a furnace for annealing. 

Placing the unannealed structure directly (optionally, first rinsing the sample) into the furnace at an 
25 elevated temperature for porosity coarsening often led to sample shape distortion or warping. The warping 
results from the inhomogeneous temperature distribution in the structure as the structure heats up, i.e., the 
structure heats up from the outside in. As a result, the porous structure undergoes coarsening first near its 
surfaces and then in its interior. The resultant surface stresses can act to distort the sample. 

Two solutions to this problem were developed and are embodiments of the present invention. The first 
30 involves the use of a rigid metal jacket cast which surrounds the structure and forces the structure to 
maintain its original shape. The cast is designed so that it is an adjustable stainless steel jacket which 
snugly surrounds the structure. A sample placed in such a jacket retained its original shape and dimensions 
during annealing at the desired temperature, in the furnace. The desired temperature ranges from 100- 
800° C. 

35 The other solution relies on the use of a temperature programmable furnace or similarly temperature 
controlled furnace. The unannealed sample is placed in the furnace at ambient temperature and the 
temperature is ramped (increased) linearly at a rate of less than 10* C min~ 1 until the desired temperature 
is attained. The desired temperature range is from 1 00-800 "C. This procedure is necessarily more time 
consuming than the procedure involving the cast where the unannealed structure is placed in the furnace at 
40 the aforementioned prescribed elevated temperature, for a prescribed time interval. 

For both procedures, the time interval for continuously maintaining the prescribed temperature ranges 
from 2 minutes to one hour and preferably 5-30 minutes. 

Both procedures can be used simultaneously to provide even less shape distortion. 

For all annealing, including in-situ annealing discussed below, the mean pore size < t> varies with time, 
45 at constant temperature, proportional to time t (0 ' 25 "°- 35) (preferably, t* 0 - 30 *) i.e. a t (, - 2W1 - 35 >. 

In one embodiment of the present invention, the annealed porous metallic structure may be reinforced 
with a metal (e.g. copper) macroscopic mesh or other design placed on the surface of the porous structure 
to provide greater structural support. The macroscopic mesh forms a composite like structure providing 
great structural integrity. 

so In another embodiment of the present invention, a step change in the porosity of the porous metallic 
structure is achieved in the following manner. Corrosion (de-alloying) is carried out to a prescribed depth 
which is some fraction of the sample thickness. Following this procedure, the sample is heat annealed 
according to the previously described procedures. The sample is removed from the furnace and once again 
subjected to corrosion for a prescribed period of time. A corrosion cell is employed which can be 

55 evacuated. The cell is evacuated after the membrane is placed back into the cell. The electrolyte is 
introduced into the cell under vacuum contacting the porous structure. This facilitates penetration of the 
electrolyte into the porous structure. A fine porous sponge forms at the interface of the coarse annealed 
porous structure and uncorroded metal. Once again the sample is placed into the furnace, but annealed at a 
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temperature lower than the first anneal. This allows the initial coarse structure to remain stable during the 
second anneal. The sample is removed and placed into the corrosion cell for further de-alloying. Again the 
annealing procedure is repeated at yet a lower temperature and a third layer of porosity is formed. The 
entire procedure may be repeated again and again thus forming a de-alloyed membrane of prescribed total 
s thickness with step changes in porosity. 

The present invention is further illustrated by the following example: 



EXAMPLE 

10 

The processes of the present invention were employed for the production of the porous metallic 
structures of the invention, for example, nano- and micro-porous gold metal membranes as follows: The 
membranes were made from a silver-gold (Ag-Au) alloy. The AE° for this system was about 0.7 Volt so 
that the de-alloying (i.e. removal of Ag) process occurs via a two-dimensional spiral process. Silver Is the 

75 more reactive element and so a minimum concentration of 60 at % Ag was required in the alloy in order to 
produce porous gold. A 1M perchloric acid electrolyte was employed owing to the high solubility of silver in 
this solution. Various concentrations of perchloric acid were tested ranging from 0.01 M to 4M. Other than 
perchloric acid, nitric acid was also tested and found to be an appropriate electrolyte. Perchloric acid was 
found to be more useful in that there was a direct correlation between the quantity of Ag removed from the 

20 alloy and the total anodic charge measured as the integral of the current-time plot The anodic current 
provided a direct measure of the number of Ag atoms removed every second. Since silver dissolves as 
Ag*, each dissolving Ag atom is responsible for a charge of ca 1.6 x 10~ 19 coul. The current was directly 
measured using the ammeter of the potentiostat. The current density is the current divided by the exposed 
area of the sample which is dissolving. A very high conductivity electrolyte was employed in order to lessen 

25 iR (voltage) drops down the microporous channels. The voltage drop greatly reduces the kinetics of the de- 
alloying process and can considerably lengthen the time required to obtain the porous solid. A 1M 
perchloric acid concentration was used because it represents a good balance between the high conductivity 
requirement and safety, i.e., perchloric acid can be an explosive mixture in concentrations between 4-1 0M 
at elevated temperatures. The employed electrochemical potential varied between 0.7 Volt and 1 .4 Volt with 

30 respect to a saturated calomel electrode (SCE). The optimal voltage vis-a-vis efficiency was ca 1.1 Volt 
(SCE), as shown in FIGURE 1. 

The Ag-Au alloys were melted in the form of 1 cm. rods or ingots under vacuum using standard 
procedures. The ingots were furnace cooled and then swaged so that the rod diameter was reduced by 
approximately 10%. Following this cold working, the alloys were placed in a furnace under vacuum and 

35 annealed at 900 *C for 120 hrs. This procedure substantially removed the dendritic microstructure initially 
present in the as cast material, which while not necessary, provided more uniform strength of the resulting 
porous structure. 

The behavior of Ag-Au alloys was examined in the range of compositions of 10-90 at % Ag and it was 
determined that the optimum composition range for gold membrane and sponge (generally, porous 

40 structure) production was 60-80 at % Ag. These results are in agreement with the previously discussed 
percolation model of de-alloying developed by the present inventors. FIGURE 1 shows polarization curves 
for selected compositions in the range examined and demonstrates the significant alteration in behavior 
which occurs near 60 at % Ag. For alloy compositions containing less than ca 60 at % Ag, the alloy 
polarization curve is similar to the polarization curve for pure gold. For alloy compositions with the Ag 

45 content greater than 60 at %, the current density rises significantly above the 0.01 mA-cm 2 level which is 
characteristic of pure gold and lower (<60 at %) silver containing alloys in the perchloric acid electrolyte. 
Note that the sustained plateau current (at ca 1.1 V (SCE)) of the alloys does not simply scale with the silver 
concentration in the alloy. This is due to the fact that the selective corrosion of silver occurs by a 
percolation like process. The size of the porosity that develops from the silver removal controls the rate of 

so the subsequent silver removal and is a function of the starting silver content in the alloy. The lower the 
initial silver content the smaller the porosity which develops in the as-corroded structure. Fine porosities will 
slow down the kinetics of mass transport in the electrolyte and increase iR (voltage) drops within the pores 
and these effects greatly reduce the kinetics of the de-alloying process. In general within the 60-80 at % Ag 
range, lower silver contents led to finer porosities in the as-corroded structure which slowed the kinetics of 

55 the de-alloying process. Near 60 at % Ag, the de-alloying process was extremely slow with sustained 
currents in the range of 0.1-1.0 mA-cm 2 . These current densities were not conducive for the production of 
thick (MOO urn) membranes. At the other extreme, the selective dissolution of the 80 at % Ag containing 
alloy resulted in gold porous structures which were extremely fragile in the as-corroded state. Scanning 
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electron microscopy showed that numerous cracks can develop in gold porous structures produced from 

alloys close to this composition. 

In general, it was found that compositions close to 60 at % Ag were best for the thin (<10 urn) fine 

porosity (10 nm) membranes and that compositions close to but less than 80 at % Ag were best for the 
5 production of thick membranes of large average porosity ( about 1 0 u,m). 

An optimal composition was found for which a large variety of very stable porous gold membranes 

could be produced. This composition is Ago.74Auo.26 and this alloy polarization curve is shown in Figure 1. 

Samples of this composition could be machined prior to the corrosion procedure into any prescribed shape 

(e.g., cylindrical, disc, etc.,) which maintained very close tolerances during the production of the membrane. 
10 Various electrode configurations were examined for mechanical stability, i.e., a corroded (at 25 " C) 

unsupported gold sponge was very susceptible to collapse under its own weight. It was determined that 

supporting the Ag-Au electrode on a large mechanically stiff very flat metal sheet (3), as shown in Figure 2, 

made of platinum or another stable metal which serves as part of the working electrode, completely 

eliminated the problem of poor mechanical stability. The electrochemical cell configuration is shown in 
75 Figure 2. Structural support of the membrane during the corrosion process was found to benefit stability. 

It was found that in-situ annealing can be accomplished by performing the de-alloying at about 50-90* C 

and preferably 90 * C. In-situ annealing increased the as-corroded ductility of the resulting porous structures. 

Moreover, in-situ annealing enabled the manufacture of membranes with geometries other than planar which 

were stable mechanically. 

20 The Ago.74Auo.26 electrode (1) was configured for de-alloying as shown in Figure 2., in a 1.0M perchloric 
acid electrolyte. The Ag 0.74AU0.26 alloy was placed on the support plate which was connected as the 
working electrode to a potentiostat. The counterelectrode (2) material in this configuration was platinum but 
. any suitable metal (e.g. stainless steel) or graphite would serve equally well. The electrochemical potential 
was set at ca 1.1 Volt (SCE) and maintained at that potential until the anodic current density dropped below 

25 the level of 10 microamps-cm~ 2 . At this current density the anode composition was >99% gold forming 
porous structure (4). This procedure was performed at temperatures in the range 1O-90*C to form 
membranes (depending upon the desired characteristics of the membrane). Following the de-alloying 
treatment, the anode, now composed of nanoporous gold, can be transferred to a temperature controlled 
furnace for the annealing (coarsening) procedure. 

30 Membranes of dimensions 10mm x 10mm x i.5mm thick were fabricated from an ingot of Ago.74Auo.26 
alloy using the ambient temperature selective corrosion process described above. Subsequent samples 
were encapsulated in a cast stainless steel jacket and inserted into the annealing furnace at a prescribed 
temperature for a time of 10 minutes. The resultant membrane characteristics are listed in the table below. 

35 
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TABLE: 



Characteristics of membranes manufactured from a Ag 0 ,74Au 0 ^B 
alioy. 

The Table below lists the resulting characteristics of the 
membranes obtained. The abbreviations in the table refer to the 
following: E, the elastic Young's modulus of the menbrane; S the 
fracture stress; < 1> the average pore size of the membrane; and 
A, the average pore area; 



Temp(°C) 


A<um 2 ) 


< 7>(um) 


E(MPa) 


S(MPa) 


25 (room temp.) 


6x1 0~ 6 


0.002* 




Mk 


100 


2.5x1 0~ 5 


0.005 


240 


0.55 


200 


0.0004 


0.02 


190 


0.65 


300 


0.002 


0.04 


160 


0.80 


400 


0.01 


0.10 


100 


1.10 


500 


0.02 


0.14 


100 


1.30 


600 


0.07 


0.26 


97 


1.45 


700 


0.20 


0.44 


55 


1.50 


800 


1.53 


| 1.25 


47 


3.00 



•Estimated by linear extrapolation 
"Not determined 



The void size distributions for the samples annealed at 800, 700, 600, 500 and 400* C are shown in 
Figure 3, and the aspect ratio distributions for these samples are shown in Figure 4. The distributions 
almost superpose on one another indicating the self-similar nature of the coarsening process operative 
during annealing. Thus the uniformity of the pore size is assured within the range of temperatures 
examined. 

The mechanical properties of the membranes are listed in the table above. In general the smaller pore 
size membranes have a larger Young's modulus and a smaller fracture stress than the larger pore size 
membranes. It was also determined that the toughness or ductility of the membrane was a function of the 
ratio of the membrane thickness, W, to the average pore size. For membranes with W/<7> less than 
10,000 the membranes were tough and ductile, while for membranes with W/</> significantly greater than 
10,000 the membranes were brittle. It was found that the very small pore size membranes can be 
toughened or ductilized by decreasing the thickness of the membranes such that the value of W/<r/> is 
less than 10,000. 

While only a few exemplary embodiments of the present invention have been described in detail, those 
skilled in the art will recognize that the invention is not restricted to said embodiments but is defined in the 
following claims. 

Claims 

1. A process for producing micro- and nano-porous metallic structures comprising contacting an alloy of 
at least two metals having different electrochemical activity in predetermined proportions with an electrolyte 
and applying an electrochemical potential to said alloy wherein said alloy acts as the anode until all or a 
portion of the more electrochemically active metal is dissolved from the alioy leaving porous metal. 

2. A process as claimed in Claim 1, wherein said alloy is a binary alloy containing 60-80% of the less 
noble metal and 40-20% of the the more noble metal. 

3- A process as claimed in Claim 1 or Claim 2, wherein said alloy is a binary alloy of silver or gold. 

4. A process as claimed in Claim 3, wherein silver forms 74% and gold forms 26% of said binary alloy. 

5. A process as claimed in Claim 3 or Claim 4, wherein said electrolyte is perchloric acid or nitric acid. 

6. A process as claimed in any one of the preceding claims, wherein said porous metal is heated for a 
predetermined period of time to anneal said porous metal. 

7. A process as claimed in Claim 6, wherein the temperature of said electrolyte is maintained at 50- 
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90 *C, for a predetermined period of time, for in-situ annealing. 

8. A process as claimed in any one of the preceding claims, wherein selected areas of said alloy are 
masked before contacting with said electrolyte, to prevent or reduce subsequent dissolution of the less 
noble metal from said covered areas, the undissolved areas providing support to the resulting porous metal. 

9. A process as claimed in any one of the preceding claims, wherein the less noble metal is dissolved 
from only the surface of said alloy to a predetermined depth. 

10. A process as claimed in Claim 9, wherein: 

(a) said porous metal is heated for a predetermined period of time, to anneal said porous metal; 

(b) the annealed porous metal is placed into a corrosion cell and said cell evacuated; 

(c) electrolyte is introduced into said evacuated cell, contacting the porous metal; 

(d) an electrochemical potential is applied to said alloy wherein said metal acts as an anode until 
additional less noble metal is dissolved from beneath the porous metal leaving fine porous metal at the 
interface of the annealed porous structure and undissolved metal; 

(e) said porous metal of step (d) is heated at a lower temperature than annealing step (a) for a 
predetermined period of time, to anneal said porous metal; and 

(f) optionally, repeating steps (a)-{e) for a predetermined ("n") number of cycles to produce a porous 
metal of prescribed thickness with "n" number of step changes in porosity with depth from the surface of 
said porous metal. 

11. A porous metal structure obtainable by a process as claimed in Claim 1 and having substantially 
uniform and interconnected porosity, with pore sizes ranging from 3 nanometers to 1 0 micrometers. 

12. A porous metal structure as claimed in Claim 11 comprising a binary alloy substantially formed of 
gold with less than 1 % silver. 

13. A porous metal structure as claimed in Claim 11 or Claim 12, wherein said structure has a thickness 
ranging from 0.0001 to 10mm. 

14. A porous metal structure as claimed in any one of Claims 11 to 13, wherein the degree of porosity 
changes with depth in a step fashion from the surface of said structure. 
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Dealloying is a common corrosion process during which an alloy 
is 'parted* by the selective dissolution of the most electro- 
chemically active of its elements. This process results in the 
formation of a n an op or o us sponge composed almost entirely of 
the more noble alloy constituents 1 . Although considerable atten- 
tion has been devoted to the morphological aspects of the deal- 
loying process, its underlying physical mechanism has remained 
unclear 2 . Here we propose a continuum model that is fully 
consistent with experiments and theoretical simulations of alloy 
dissolution, and demonstrate that nanoporosity in metals is due 
to an intrinsic dynamical pattern formation process. That is, 
pores form because the more noble atoms are chemically driven 
to aggregate into two-dimensional clusters by a phase separation 
process (spinodal decomposition) at the solid-electrolyte inter- 
face, and the surface area continuously increases owing to etching. 
Together, these processes evolve porosity with a characteristic 
length scale predicted by our continuum model. We expect that 
chemically tailored nanoporous gold made by dealloying Ag- Au 



should be suitable for sensor applications, particularly in a 
biomaterials context 

Selective dissolution has a long history 5 . For example, the 
chemical treatment known as depletion gilding selectively dissolves 
a non-gold element near the surface of a less-expensive alloy such as 
Au-Cu, leaving behind a surface of pure gold. Early Andean gold- 
smiths used this technique to enhance the surfaces of their artefacts 4 . 
In this century, selective dissolution has been primarily examined in 
the context of corrosion. It is observed in technologically important 
alloy systems, notably brasses, stainless steels, and Cu-Al alloys 1 * 6 . 
The mechanical properties of a porous overlayer are very different 
from the bulk alloy on which it sits, leading to brittle crack 
propagation, stress corrosion cracking, and other undesirable 
materials failure 7 . Figure 1 shows the prototypical dealloyed micro- 
structure, that of nanoporous gold. Early notions considered 
porosity as a hidden microstructure revealed by etching, but 
diffraction experiments showed that no pre-existing length scale 
exists before acid attack on single-phase alloys* 1 *. Later ideas 
considered the influence of percolating clusters within the solid 
solution of the alloy, but models railed to yield behaviour consistent 
with experiment ,0 ' n . 

The following argument illustrates the fundamental obstacle to 
understanding porosity formation during dealloying: consider a 
silver-gold alloy in an electrolyte under conditions where silver 
dissolves and gold is inert. Initially, silver will be dissolved from 
surface sites such as terraces or steps. Gold atoms should accumu- 
late on the surface and locally block further dissolution. For an alloy 
containing 10% gold, it might be expected that dissolution would 
stop or be significantly retarded after about 10 monolayers of the 
alloy have been dissolved. 

A complete model of selective dissolution needs to be multi-scale, 
involving the kinetics of dissolution, surface diffusion, and mass 
transport through the bulk of both alloy and electrolyte. Because 




Figure 1 Scanning electron micrographs of nanoporous gold made by selective 
dissolution of siter from Ag-Au alloys immersed In nitric add under free corrosion 
conditions, a, Cross-section of deaBoyed Au^^Ag^ (atom%) thin film. b. Plan view ot 
deaiioyed Au M *A9 74% (atom%). The porosity is open, and the ligament spacings shown in 
b are of the order oi lOnrn; spacings as small as 5 nm have been observed. 
Measurements oi the surface area of nanoporous gold are of the order of order 2nrV* 
(refs 24 , 25), comparable to commercial supported catalysts. 
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mass transport through the bulk of the growing phase (the electro- 
Ivte) is always a stabilizing influence", and because mass transport 
through the bulk of the dissolving phase appears too . slow to be 
significant, we hypothesised that the morphology-determining 
physical process is confined to the interface region between the 
alloy and the electrolyte. To test this, wc developed a kinetic .Monte 
Carlo model to simulate Ag-Au dealloying. including only difosion 
of silver and gold and dissolution of silver". We found that this 
model was able to reproduce all relevant experimental trends 
characteristic of dealloying, both morphological and kinetic 

Figure 2 shows a simulated porous structure with hgament widths 
of 2-5 nm. Our simulations were successful in modelling the 
nanoporous morphology, and also in modelling ^J^ c 
behaviour of the dissolution current versus overpotendaL It is a 
well characterized feature of alloy dissolution that as the over- 
potential is increased (usually at rates of the order of a few milhvo ts 
per second), the dissolution current of ions from the alloy stays at a 
low levd until a bulk-compos^^ 

is reached, at which point this current rises rapidly 1 Figure 3 
shows simulated and experimental polarization curves for different 
alloy compositions. There is dear observation of a composition- 
dependent V c . To our knowledge, this is the first simulation model 
to produce such behaviour, suggesting that we have found a 
minimum set of physical processes to include in any model for 




Rgure 2 Simulated nanoporous gold. The simulation model was as fellows: a bond- 
breaking model was used for diffusion; atoms with N near neighbours diffused wtm rate 



^ ppexphNe/feT), where ( is a bond energy and vo = 10 13 s" 



. Dissolution rates were 



consistent with the Butier-Voimer (BV) equation In the hlghMJrhrlng-force Tatel regime; 
the dissolution rate fe, for a silver atom with N near neighbours was written 35*u/= 
v expH* - 0 W). wnere v £ = 1 0 4 s- 1 is an attempt frequency determined by the 
exchange-current density m the BV equation and * Is the overpotential. For the figure, 
$ st 1.75eV,«/AeT"*5.51. 



alloy dissolution. . c 

The simulations reveal the following qualitative picture of 
porosity formation. The process starts with the dissolution of a 
single silver atom on a flat alloy surface of close-packed (111) 
orientation, leaving behind a temce vacancy. The atoms coordinat- 
ing this vacancy have fewer lateral near-neighbours than other silver 
atoms in the tenace, and are thus more susceptible to dissolution. 
As a result, the entire tcnace is 'stripped', leaving behind gold atoms 
with no lateral coordination Oadatoms*). Before the next layer is 
attacked, these gold adatoms diffuse about and start to agglomerate 
into islands. As a result, rather than a uniform diffuse layer of gold 
spread over the surface, the surface is comprised of two distinct 
kinds of regions— pure gold dusters that locally passrvate the 
surface, and patches of un-dealloyed material exposed to electrolyte 
When silver atoms in these patches dissolve, more gold adatoms are 
released onto the surface. These adatoms diffuse to the gold dusters 
left over from dissolution of previous layers, continuing to leave un- 
dealloyed material exposed to electrolyte. In the early stages, these 
gold dusters arc mounds that are gold-rich at their peaks and that 
have alloy composition at their bases. These mounds get undercut, 
increasing the surface area that gold must cover to bnng about 
passivation. Ultimately, this leads to pit formation and porosity. 

Central to this description is the the coalescence of gold adatoms 
into stable dusters. The spacing between these 'islands' in the initial 
stages of dissolution is dose to the spacing between ligaments in the 
final porous structure. The physical reason for this coalescence can 
be understood by considering the gold adatoms to be one compo- 
nent of a two-component solution of gold and 'electrolyte' confined 
to the monolayer-thick interracial layer sitting on top of un- 
dealloyed material. We modelled the thermodynamics of the inter 
facial layer as a regular solution 1 *, and found the solubility of gold in 
electrolyte within the interfacial layer to be of the order of 1CT per 
site (see Methods). This solubility 'may be interpreted as the 
'equilibrium concentration of gold anatoms* on the surface of the 
alloy— in the absence of etching, it represents a dynamic equilib- 
rium of adatoms resulting from their two-dimensional evaporation 
from step edges onto terraces and their subsequent rccondensation. 

In contrast to the equilibrium condition, rapidly stripping a 
terrace of sflver atoms leaves gold adatoms with a local site 
occupancy fraction equal to that in the bulk, typically 10-40%— 
far above their equilibrated concentration of 10 per site. Thus, 
there is an extremely strong driving force for gold adatoms to 
condense onto nearby gold-rich dusters. In feet, regions of surface 
with high enough supcrsaturation of gold adatoms sit within the 
spinodal 1 , a special segment of the curve representing free energy / 
of a spatially uniform layer versus gold concentration C for which 
d 7 f fdC? < 0. Within the spinodal, composition fluctuations of 
iiuinitesimal amplitude lead to a lower overall free energy for the 
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Rgure 3 Comparison of experimental and simulated current-potential behaviour 
a Currem-polential behaviour lor varying Ag-Au alloy expositions (atom* Au) 



current-potential behaviour ot Ag-Au alloys, c, Ccrnparison ol experimental dine) and 
simulated (triangles) critical potentials; the zero ot overpotential has been set equal to the 
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system, 'and involve atomic dififusion against concentration gradi- 
ents (the 'uphill diffusion* process through which gold condenses 
onto nearby clusters), that is, the system is inherently unstable and 
wul spontaneously phase-separate. But fluctuations of long length 
scale grow slowly due to the required diffusion times, and short- 
length-scale fluctuations create much energetically unfavourable 
incipient interface between the phases, inhibiting their growth. 
Hence, phase separation is manifested most rapidly at an inter- 
mediate length scale that roughly corresponds to the spacing 
between the observed gold-rich clusters. This effect is known as 
spmodal decomposition 1 "'. As porosity forms, the decomposition 
is occurring on a non-flat, non-uniform surface with continuously 
increasing surface area. 

In our model, the motion of the alloy-electrolyte interface is fully 
described mathematocally by the flux of diffusing adatoms the 
velocity of the lnt erface normal to itself v„, and the concentration 
accumulauon rate dCJdt, aD of which are interrelated and vary with 
p0S ^°" curve of the interface (for detailed derivations, 

set Methods). For J„ we used a model for diffusion during spinodal 
decomposmon known as the Cahti-Hilliard equation' 7 The 
normal velocity depends on C and also on the local curvature k 
trough capdlary effects". The time evolution of C is uniquely 
determined by the local mass-conservation condition 



aadt=v B C„-v.KC-V-/ s (i) 

where Q is :the bulk gold concentration. This condition is analogous 

L^I^T^ °^ ° r S0lute 'M*"** in boundary- 
i of ,f ^«non», with two important different 
U ) the interfaaal layer thickness is constant along the interface, and 
is microscopic, rather than being a spatially varying macroscopic 

tT^Z ^'^ ^ SUrfeC£ a ^ a »°nP-L inhere^ 
Ae Cahn-HiU,ard form for J s is essential for porosity formation. 
Simple ('downhill') surface diffusion (/ s = - Ds VQ yields a^ 




Figure 4 Simulated evolution of an artificial pit in Au )0% A 990% (atom*), * = 1 8 eV 
cross-sectons along the(111) plane defined by the yellow line In a are shown below each 
Ptan WB w. Tha UU condition is a surface fufly parted with goSS^wZt 

«* P.I e, After 1 0s. a pold cluster has nucleated in me centre of the pft. „, M 100s £ 
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innially unstable mterface that passivates quickly, before well, 
formed pores have a chance to develop. Cahn-Hilliard diffusion 
also dominates capillarity-driven surface diffusion-the effect 
usually mcorporated into interface evolution equations" 
We performed numerical integration of equation (1) usi„ E , 

matched those used u, the kinetic Monte Carlo simulations. Wc 
observed, as expected, the evolution of gold clusters separated bra 
charactenstic spacing \. An analytic expression for \ can be found 
by a time-dependent linear stability analysis of equation (1) that 
takes .mo account the slow increase of gold concentration into the 

a, dt ^ cause . hc , s P a ^ [P^od with the largest amplification 
rate depends sensitively on the gold concentration. Specifically, this 
spatial penod decreases sharply as the concentration increases past a 
direshold concentration for instability that corresponds to the 
Pinodal point aJ/dC* = 0, and the interface is stable for concert 
S 0w , thls ^^old. We find a maximally unstable spatial 
penod that scales as X « (A/V,) 1 ", where V, is the velocity of a flat 
afloy surface wth no gold accumuJatcd j t ^ . . 

quahtanve agreement with both kinetic Monte Cario simulations 
and experiments, both of which show that the characteristic length 

lho^ PC T? ^ Wi * increasin 8 "rising force. A more 
elaborate analysis mcorporating nonlinear effects, however, remains 

needed for a detailed quantitative comparison 

There is an analog,- between this result, applicable to etching, and 

^o-dimensional island nudeation during submonolayer vapour 

dSSJ eP 0tt Namd> :' fa * e "fr sta «« ° f etching, the 
dissoluOon process is analogous to deposition of gold; in both 
processes, adatoms are added to the surface where they are fr«£ 

bet Tn,^ 1ShndS - "* ° f de P«»ition has 

been studied using rate equations that describe an aggregation 
process where adatoms stick together irreversibly, fa th&luL 
it is a well-known result that the island.spadng scales as UVf)*, 
where the deposition rate F is the direct analogue of the surftce 
velocity m etchmg, and the Exponent ^depends on details of the 
aggregate process (see ref. 22 and references therein). That these 
results are limned to irreversible aggregation during deposition and 
our analysis is for reversible aggregation during etching suggests die 
«,stence of umversal scaling laws for aggregation* lhafdo not 

JxoSeT TeVenM>tY ° f * Uck * ereof in "PP^e 
Our kinetic Monte Carlo simulation eluddates the later stages of 
morphological evolution, and the mechanism by which three- 
dimensional porosity evolves. We highlight the features of this 
process by showing in Fig. 4 a simulation of an artificial pit in» 
oAerwue folly passrvated surface. When the pit reaches £ 
dept£ its surface area has increased sufficiendy that a new gold 
duster nucleates When this happens, the pore splits into multiple 

S V aCh W,th 3 SmaUer surfece area *» ' * P^ent. Th«e 
child pns continue to penetrate into the bulk, increasing their 
surface area, nudeating new dusters, spawning new pits, and so on 

• Zt % f' ^"• d j IDCnSi0nal P° rous »™*~ evolvU suTa 
those illustrated m Figs 1 and 2. g 

Methods 

.eperin. m«,,ck» en. rgi e, to^JZSZX^^fc^ 

»^ n 7»''-o^.nrfo» lhe ,. g .„ c.oflhcinUSh^t 

£bZ Z IT ?? ^ ,b f" U,C ««"»«««»• «» «i">«I«»„ cLdiU.^. 
t^*. - tV » -,. lht .imubiion bond energy,, deicribed i„ F j.. 2 1. r= 600 K 

n«P to doubled form* and , minimum al c „ , 0 -' " ri « ' 

tT? ? , ?.',r ,,U r i S 0fg0,d eleC,rol »" e <" d Tte "»)• 
ihe i^hn-H,ll u ,J J.f lu , ion equation ij/, = -MtCM^/Z.^^C 4 JM(CVP»C. 
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Here 453M(Q U a mobility, w U ihe so-called gradient energy coefficient, and the 
oradienU are laken with respect to arc length. Tlie ft r* term on the right-hand side 
describe the chemical effect leading to phase separation within the spinodal; the second 
1 term describes the effect that damps short -wavelength fluctuations. The mobdily « 
proportional to the surface diffusMty D* and is given by M(C) = (IVWW1 -c). The 
mobility is peaked for c « 0.5, and is *uo for c = 0 and t = I (atoms do not diffuse in pure 
phases because there are no vacancies in our model). The normal velocity is given by 
\ (O ^ V(QU - iyCm,TU), where -y is the surface free energy and V(0 « «" 
interface response function, equal to the velocity of a nat surface covered wtlh a 
concentration C of gold. We find in both simulation and experiment that the interlace 
Z?kL weB by the functional form V(Q * W)exp(-aC-). where * is the 
overpotential and C b a constant. E^erimenutty, the gold accumulation can be inferred 
by integrating the dissolution current versus lime at faed overpotenual; « is necessary to 
use an overpotenual that « low enough to ensure that the surface remains planar (that u, 
porosity does not form) and also to catch the short initial transient rise in c^entaxiQver 
atoms art pulled from the first few monolayers. This particular form for the Interface 
response function is quite curious. Naively, one might expect that the local interface 
velocity would be proportional to the local concentration of .flver e^ed to the 
electrolyte that u, V(C) = (1 -t). However, the decaying exponential form suggests that 
there is an evolving distribution of holes opening and dosing within the intertacia! region. 
coniroBuu the accumulation rate- 

Physical, thx mass conservation condition (equation (1)) is the statement tot the 
total number OAs of gold atoms in a length As of interface with Uteral width b can 
change as a result of three distinct effects that correspond to the three terms on the 
rUhtW-a.de of equation (l)r the accumulation of gold atoms into the > mtafeaal 
liver from the solid being dissolved; the local stretching of the interface (3As/ar * V 
A which can either increase or decrease C depending on whether the sobd is conave 
(« > 0) or convex (a < 0); and the motion of atoms along the interface driven by the 
surface diffusion flux Jy 
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It has been proposed that during Pleistocene glaciations, an ice 
cap of 1 kilometre or greater thickness covered the Arctic 
Ocean 1 - 1 . This notion contrasts with the prevailing view that the 
Arctic Ocean was covered only by perennial sea ice with scattered 
icebergs^ 6 . Detailed mapping of the ocean floor is the best means 
to resolve this issue. Although sea-floor imagery has been used to 
reconstruct the glacial history of the Antarctic shelf , little data 
have been collected in the Arctic Ocean because of operational 
constraints 10 -". The use of a geophysical mapping system during 
the submarine SC1CEX expediUon in 1999 u provided the oppor- 
tunity to perform such an investigation over a large portion of the 
Arctic Ocean. Here we analyse backscatter images and sub-bottom 
profiler records obtained during this expedition from depths as 
great as 1 kilometre. These records show multiple bedforms 
indicative of glacial scouring and moulding of sea floor, combined 
with large-scale erosion of submarine ridge crests. These distinct 
glaciogenic features demonstrate that immense, Antarctic-type 
ice shelves up to 1 kilometre thick and hundreds of kilometres 
long existed in the Arctic Ocean during Pleistocene glaciations. 

The central Arctic Ocean contains relatively shallow areas (water 
depths < 1,000 rn; see Fig. 1) on Yermak plateau, Lomonosov ridge 
and Chukchi borderland— which includes Chukchi plateau, Chuk- 
chi rise and Northwind ridge. During the SQCEX-99 expedition, 
conducted on the nuclear-powered submarine USS Hawkbill, 
shallow sea-floor areas were targeted for mapping to detect gjado- 
genie bedforms. Sea-floor images (collected using a submarine- 
mounted 12-kHz swath bathymetry and sidescan sonar 0 } from the 
Chukchi borderland and the Lomonosov ridge show a variety of 
bedforms, including random or subparallel scours, parallel linea- 
tions, and transverse ridges. On the records from the chirp sub- 
bottom profiler, these bedforms are associated with planed ridge 
crests with rough rnicrorclief and obvious angular unconformities 
cut into the stratified sediments. 

Randomly oriented furrows, typically <100-m wide and up to 
30-m deep, densely cover the shallowest, <400-m-decp portions of 
sea floor on the Chukchi borderland and adjacent continental 
margin (Fig. 2a). Isolated larger scours up to 700-m wide and 
over 10-km long occur as deep as 500 m, Even greater depths, 
exceeding 900 m, are attained by closely spaced, subparallel scours 
on the Lomonosov ridge. Sea-floor scours are known to be formed 
by the drift of icebergs and pack-ice ridges". At present, icebergs in 
the Arctic Ocean have at most 50-m draughts 14 , whereas icebergs off 
Antarctica and Greenland reach depths of 500-550 m (refs 15, 16). 
The largest depths of gouged sea floor, extending to 850 m, have 
been reported from the Yermak plateau' 0 , matching the depth of 
sours on the Lomonosov ridge. 
Below the depth range of dense scouring, the sea floor exhibits 
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